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ABSTRACT 

We investigate the impact of instrumental systematic errors on the potential of cos- 
mic microwave background polarization experiments targeting primordial _B-modes. 
To do so, we introduce spin- weighted Muller matrix-valued fields describing the linear 
response of the imperfect optical system and receiver, and give a careful discussion 
of the behaviour of the induced systematic effects under rotation of the instrument. 
We give the correspondence between the matrix components and known optical and 
receiver imperfections, and compare the likely performance of pseudo-correlation re- 
ceivers and those that modulate the polarization with a half-wave plate. The latter is 
shown to have the significant advantage of not coupling the total intensity into polar- 
ization for perfect optics, but potential effects like optical distortions that may be in- 
troduced by the quasi-optical wave plate warrant further investigation. A fast method 
for tolerancing time-invariant systematic effects is presented, which propagates errors 
through to power spectra and cosmological parameters. The method extends previ- 
ous studies to an arbitrary scan strategy, and eliminates the need for time-consuming 
Monte-Carlo simulations in the early phases of instrument and survey design. We il- 
lustrate the method with both simple parametrized forms for the systematics and with 
beams based on physical-optics simulations. Example results are given in the context 
of next-generation experiments targeting tensor-to-scalar ratios r ~ 0.01. 

Key words: cosmic microwave background - methods: analytical - methods: numer- 
ical. 



1 INTRODUCTION 

There is currently a great deal of interest and activity in the 
field of cosmic microwave background (CMB) polarimetry. 
Several recent experiments (Kovac et al. 2002; Readhead 
et al. 2004; Barkats et al. 2005; Montroy et al. 2006; Page 
et al. 2006; Wu et al. 2007) have reported detections of CMB 
polarization at a level consistent with predictions in simple 
cold dark matter (plus A) models fit to the temperature 
anisotropy data. In standard adiabatic models, polarization 
measurements promise a further tightening of constraints 
inferred from the temperature anisotropies, and the break- 
ing of important degeneracies (Zaldarriaga, Spergel & Sel- 
jak 1997). In particular, with the recent release of the three- 
year Wilkinson Microwave Anisotropy Probe ( WMAP) data, 
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the optical depth to reionization has been constrained to 
r = 0.09 ± 0.03 (Spergel et al. 2006) with the large-angle 
polarization signal (Zaldarriaga 1997). Polarization informa- 
tion is particularly valuable in non-standard models, for ex- 
ample it can significantly improve constraints in the presence 
of isocurvature modes (Bucher, Moodley & Turok 2001), and 
be used to constrain parity-violating physics (e.g. Scanna- 
pieco & Ferreira 1997; Lue, Wang & Kamionkowski 1999; 
Feng et al. 2006) . Secondary effects, most notably weak grav- 
itational lensing (Zaldarriaga & Seljak 1998; Lewis & Challi- 
nor 2006), further encode information on the low-redshift 
universe in CMB polarization. 

Perhaps the most exciting aspect of CMB polarime- 
try is the window it may open on primordial gravitational 
waves (Kamionkowski, Kosowsky & Stebbins 1997; Zaldar- 
riaga & Seljak 1997). The curl mode (or B-mode) of polar- 
ization is not produced by linear density perturbations and 
so provides an observable signature of gravitational waves at 
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last scattering that is not limited by cosmic variance from 
the dominant density perturbations. However, the situation 
is clouded by several issues. First, the amplitude of grav- 
itational waves expected from inflation is unknown as the 
energy scale at which inflation may have occurred is not 
(yet) determined by theory. What is known is that any B- 
mode imprint will be very small: the current limit on the 
amplitude of the gravitational wave power spectrum (ex- 
pressed as a fraction of that for the density perturbation) 
is r < 0.28 at 95% confidence (Spergel et al. 2006), from 
a combination of the temperature anisotropies and galaxy- 
clustering data. This limit on r translates to an r.m.s. B- 
mode signal < 200 nK. Furthermore, gravitational lensing 
does produce B-mode polarization at second order in the 
density perturbations, with a spectrum that is almost white 
for multipoles I < 300. The amplitude ^/C ; B lons « 5/iK- 
arcmin means the lens-induced B-modes dominate the pri- 
mordial ones for r < 0.0f except on very large scales where 
the latter is enhanced by reionization. For noise levels much 
better than 5 /iK-arcmin, it will be worthwhile to try and 
clean out the large-angle B-modes of lensing by making use 
of non-Gaussianity (e.g. Hu & Okamoto 2002). Finally, po- 
larized emission from Galactic foregrounds will likely over- 
whelm the primordial B-mode signal over a large fraction 
of the sky (Page et al. 2006), and our ability to constrain 
gravitational waves will almost certainly be limited by the 
accuracy at which Galactic foregrounds can be subtracted. 

The intrinsic weakness of the polarization signal 
presents a major experimental challenge. The gradient (or 
B-mode) polarization is now determined to be at least an or- 
der of magnitude smaller than the temperature anisotropies 
for / < 1000, and the B-modes (including lensing) are in- 
directly constrained to be at least an order of magnitude 
smaller still. As well as the raw sensitivity requirements that 
this implies, any systematic errors present in B-mode instru- 
ments will have to be controlled to an unprecedented level 
of accuracy in order to ensure the signal is not fatally con- 
taminated. With this motivation, in this paper we develop 
a general framework to describe and assess the impact of 
instrumental systematic errors on experiments targeted at 
CMB polarization. 

There is a growing literature discussing the impact of 
systematic effects in CMB polarimetry; for discussions of the 
real-world issues encountered in recent surveys, see Barkats 
et al. (2005); Masi et al. (2006); Page et al. (2006); Jones 
et al. (2006). There are typically many sources of poten- 
tial systematic error. Here we shall be concerned with the 
broad class of errors that produce a time-independent resid- 
ual signal for a given instrument pointing. These include 
imperfections in the receiver and optics, but exclude impor- 
tant effects such as low-frequency noise from thermal fluctu- 
ations, readout electronics or atmospheric fluctuations, and 
pointing jitter. Low-frequency noise is best dealt with by a 
combination of signal modulation (active or by scanning), 
a well cross-linked survey, and removal during the map- 
making stage. Our aim is twofold: to give a careful anal- 
ysis of the transformation properties of various systematic 
effects under rotation of the instrument - a useful strat- 
egy for mitigating some systematics; and to provide a fast, 
semi-analytic method for tolerancing systematic effects that 
is flexible enough to deal with arbitrary scan strategies and 
removes the need for time-consuming simulations during the 



early stages of instrument and survey design. Our approach 
is similar to that of Hu, Hedman & Zaldarriaga (2003), but 
we extend their analysis in several ways. We describe the re- 
ceiver systematics with Miiller matrices and the optics with 
Miiller matrix- valued fields which allows us to deal with arbi- 
trary beam effects. To simplify the discussion of the rotation 
properties of systematic effects, we introduce spin-weighted 
Miiller matrices and perform a further decomposition of the 
beam matrix-valued fields into irreducible components. The 
latter singles out those features of the beams that cannot be 
overcome with instrument rotation, and easily reproduces 
the results of Carretti et al. (2004) for axisymmetric sys- 
tems. In addition, our analysis of the scientific impact of 
systematic effects works with an arbitrary scan strategy, and 
propagates the effects through to biases in cosmological pa- 
rameters and the increase in their random errors. Carrying 
biases through to parameters is potentially important: sys- 
tematic errors often must be suppressed further than the 
cosmic- variance limit in the CMB power spectra if they are 
to have negligible impact on the science (Efstathiou & Bond 
1999). Also, properly modelling the interaction between sys- 
tematic effects and the scan strategy is important as its ef- 
fects cannot easily be predicted, or constrained, from simpler 
scans (such as a raster scan with no implied beam rotation) . 
For example, a non-trivial scan can cause leakage from the 
dominant B-modes into B-modes that is not present for a 
raster scan, as a result of the non-local nature of the E- and 
B-mode decomposition. 

This paper is organized as follows. Section 2 introduces 
our notation and some important assumptions. Propaga- 
tion through the receiver is described in Section 3 in terms 
of spin- weighted Miiller matrices, while the optics are de- 
scribed in terms of matrix-valued fields in Section 4. We 
give illustrations for common receiver types and optical im- 
perfections, and compare their relative merits. In Section 5 
we present our method for propagating errors to biases in 
power spectra and the enhancement of their variances, and 
follow these through to cosmological parameters (specifically 
r) in Section 6. Finally, in Section 7 we apply our methods 
to set tolerances on parametrized errors and compare these 
for different scan strategies. We also examine the impact 
of realistic, simulated beam profiles. Two appendices give 
further properties of the complex Miiller matrices and of 
the beam expansion in irreducible components introduced 
in Section 4. 



2 NOTATION 

For surveys covering a small fraction of the sky we can work 
in the flat-sky limit where the CMB fields are considered on 
the tangent plane to the celestial sphere. We adopt Cartesian 
coordinates on this plane such that x points north-south 
and y points west-east; the right-handed z-direction is then 
along the line of sight. We define Stokes parameters using 
the x and —y directions for which the radiation propagation 
direction completes a right-handed triad. Our polarization 
conventions are then consistent with the IAU standards. If 
the electric field of the incident radiation at frequency u> is 
Jt^Be -1 "*), with these conventions we have 
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(EiE*) 



T + Q 
-U + iV 



-u 

T - 



-iV 

Q 



(i) 



We shall only consider quasi-monochromatic systems here 
but our discussion will still hold for broad-band systems if 
the non-ideal instrument responds in the same way to all 
in-band radiation. Clearly this ignores an important class of 
systematic effects, for example bandpass mismatch (Jarosik 
et al. 2006). 

The complex polarization P(x) = (Q + iU)(x) is spin 
—2 in the sense that under the transformation x + iy i— ► (x + 
iy)e % ^ ', which rotates the :r-axis by ip towards the negative 
j/-axis, P I— » Pe -2 '*. Decomposing P into its electric (E) 
and magnetic (B) parts, we have, in Fourier space, 



(Q±iU)(x) = 



^[E(l) T iB(l)]< 



^p2id>i il-x 



(2) 



where I = Z(cos <f>i, sin 4>i). For the cosmological examples 
we give in this paper we only consider the CMB fields on 
the sky. For observations in the quiet Galactic regions that 
will be the targets for future ground-based B-mode experi- 
ments, the CMB temperature fluctuations and i?-mode po- 
larization should dominate Galactic emission. Ignoring the 
latter should therefore be harmless for those most troubling 
systematic effects that couple temperature and E-mode po- 
larization into B. 



3 RECEIVER MULLER MATRICES 

The Miiller matrices describe the propagation of the Stokes 
parameters through the receiver element of a given observing 
system. The optical coupling to the fields on the sky requires 
a description in terms of Miiller matrix- valued fields which 
we describe in Section 4. 

Gathering the Stokes parameters in a Stokes vector s = 
(T,Q,U,V) T , we have for the observed Stokes vector 



Sobs — 
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(3) 



We adopt the convention here that Miiller matrices are al- 
ways expressed in the instrument basis. This coincides with 
the Cartesian sky basis (x and — y) when the instrument is 
in its fiducial orientation. It is convenient to work with the 
complex Miiller matrix whose elements have definite spin, 
i.e. 



fobs 
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(4) 



Here, the complex Stokes vector is p — (T, P, P* ,V) T . We 
denote the matrix on the right of equation (4) by M. Its 
components are related to those in equation (3) as follows: 
for the total intensity 



Mi 



\{M T q - iMtu), M T p* = \{M T q + iM TU ), (5) 



with analogous results for Mvp and Mvp*\ and for the po- 
larization 



and 

Mpp = 
Mpp* 



■ M uv ) + \i(M UQ - M QU ) 
- M uv ) + \i(M VQ + M QU ) 



(7) 



The components for the spin-2 polarization P* are related 
to those for P: 



Mp* T = M* PT , 
Mp* P = Mpp*, 



M P *v = M*p V , 
M p *p* = Mpp. 



(8) 



If we kept the instrument in a fixed orientation, but 
transformed to a rotated basis in describing the polariza- 
tion fields on the sky, the Miiller matrix elements would 
transform like the complex conjugate of the field appearing 
in the second index, e.g. Mpp Mppe 2 ^ so that MppP 
remained constant. If we further rotated the observed polar- 
ization from the instrument basis to the rotated sky basis, 
we would pick up an additional factor e ls ^ where s is the 
spin of the field associated with the first index. More rele- 
vant for our purposes are the transformation properties of 
the Miiller matrix under (active) rotations of the instru- 
ment. Let us rotate the instrument by tp taking x towards 
the negative j/-axis and simultaneously back-rotate the ob- 
served polarization so we are describing the measured polar- 
ization in the original sky basis. In its basis, the instrument 
sees incoming radiation with complex Stokes vector A(i/>)p, 
where A(ip) = diag(l,e~ l *,e 1), so the observed polar- 
ization on the sky basis is p obs (V0 = A'(^)MA(i|i)p. Note 
that only the diagonal elements of M are invariant under 
M i-> A t (V')MA(i/'). 

In the case of an ideal instrument, M is equal to the 
identity matrix and any systematic errors that affect the 
Stokes parameters will lead to small perturbations from this. 
It is convenient to introduce a more concise notation for 
those elements describing coupling to P and P* (Hu et al. 
2003): M PT = 71 + «72, M PP = 1 + a + 2iu>, M PP * = 
fi + if2 and Mpv = wi + iw2, where all parameters on the 
right-hand sides are real. For the perturbation to the linear 
polarization, we then have 

5(Q±iU) = {a±2ico)(Q±iU) + {f 1 ±if 2 )(QTiU) 

+ (71 ±iyz)T+ (wi ±iw 2 )V. (9) 

These parameters are defined in the fiducial basis, and de- 
scribe a miscalibration of the polarization amplitude, a, a 
rotation of the polarization orientation, to, transformations 
between the two polarization spin states, fi and f 2 , leakage 
from total intensity to Q with 71 and to U with 72, and 
leakage from circular polarization with mi and w 2 . The cir- 
cular leakage terms were not considered in Hu et al. (2003) . 
Under a rotation of the instrument as described above, we 
find 



S(Q±iU) = 



(a±2iu))(Q±iU) 
+ (fi±if 2 )(QTiU)e 
+ (71 ± n 2 )Te ±2ii> 
+ (wi ±tW2)Ve ±2z ^ 



±4iip 



(10) 



Mp T = Mq T + iMur, Mp V = Mq V + iMuv, 



(6) 



This suggests that the /, 7 and w errors can be controlled 
with instrument rotation as these terms have different spin 
properties to the fields they are perturbing. For example, if 
we average observations of a point over all possible orienta- 
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Figure 1. Block diagrams for a pseudo-correlation receiver (left) 
and a rotating half- wave-plate receiver (right). For the pseudo- 
correlation receiver, the incident radiation is split into two or- 
thogonal components by the ortho-mode transducer (OMT), and 
then propagates through a 90° hybridizer. A time dependent 
phase shift is introduced along one arm, and the radiation passes 
through a second 90° hybridizer before being detected. For the 
rotating half-wave-plate receiver, the incident radiation simply 
passes through a rotating half-wave plate before being split by 
the OMT and detected. 



tions these terms disappear. We look at controlling system- 
atics with beam rotation in more detail in Section 4.1. 



3.1 Receiver errors 

In this section we give examples of the polarization Miiller 
matrix elements for two common polarimeters that have 
particular relevance to CMB polarimetry: the pseudo- 
correlation receiver and the rotating half-wave-plate re- 
ceiver. Block diagrams for these receivers are shown in Fig. 
1. 

Quite generally, the propagation of radiation through a 
receiver can be described by a Jones matrix, J, such that 
the electric field after passing through the receiver, £/ rcc , is 



SE, 



(ii) 



where E is the incident electric field. In this section only, the 
elements of E are the complex amplitudes of the two linear 
polarizations, A and B, which, for an ideal optical system, 
couple in the far-field to the x and —y components of the 
electric field of the incident radiation; see Section 4. For a 
receiver with several components, the Jones matrix of the 
receiver is the product of the matrices for each component 
provided that reflections can be ignored. The Miiller matrix 
for the receiver can be found from the relations s obs = Ms 



T + Q U + iV 
U-iV T-Q 



.(T + Q U + iV \ t 

' U-iV T-Q J ■ 

(12) 



For the pseudo-correlation receiver, the ideal Jones ma- 
trix is 

Jpc = Jhybrid,2 Jphasc Jhybrid.l Jomt 
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71 



V2\ 



i 
1 



1 - e- 
i(l + e ivt ) 





1 
1 

_! +e iv>t 



(13) 



where ipt is the time-dependent phase shift, assumed to be 
continuous here. After passing through the receiver, detec- 
tors measure the power in the two components, E TCCi a and 
£rcc,s and, in the ideal, noiseless case, their outputs for a 
particular pointing are 

Di = (\E ICC:A \' 2 } = \{T -Qcostpt- Usaupt) 

=- <|Srcc,s| 2 ) = \{T + Q cos ipt + U sin <pt). (14) 



D 



The total intensity and linear polarization Stokes parame- 
ters can be recovered by taking the sum and difference of 
the detector outputs and demodulating. If we assume that 
the modulation frequency is much higher than the maxi- 
mum frequency in the incident signal (which is determined 
by the scan speed and resolution), we can approximate the 
demodulation step by 



Qobs 



U 



1 

A 



-r- / (D 1 +D 2 )dt 



/ (D2 — D\) cos (fit dt 
Jo 

/ (D 2 - Di) sin <ptdt, 
Jo 



(15) 



where A is long compared to l/<p. We can now follow system- 
atic errors introduced in the Jones matrices through to the 
observed Stokes parameters. We parametrize the systematic 
errors in the various receiver components as follows: 



l+3i 
tie z 



ere" 1 
(1+32K Q 



Jphasc 



J hybrid, j 



1 

7! 





e i(ipt+S<j>) 
(1 + ^) 

*(i + c 3 ; 



i(l + Bi) 

(i + ^O 



(16) 



where j = 1, 2 labels the two hybridizers. Each parameter 
corresponds to a potential, physical systematic error, for ex- 
ample, gi and <?2 represent gain errors in the two arms of 
the OMT. Following the same process as in the ideal case, 
we find the effect of these errors on the observed Stokes pa- 
rameters. Assuming that the systematic errors do not vary 
significantly over the time of the observation, and expand- 
ing to linear order, the errors are related to the parameters 
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introduced in equation (9) by 

a = \(Ai +Bi + Ci + Di + A 2 + B 2 + C 2 + D 2 ) 
+9i + 92 

2lj = i(oi + bi - ci - di + a 2 + c 2 - 62 - d 2 ) 
+e 2 cos #2 — £i cos 61 — 5</> 

71 = 3i -32 + 5(^1 + C*i - Si - Di) 

72 = ei cos #1+62 cos 6 2 + 5 (ai + di — 61 — ci ) 
/1 = 

/2 = 

wi = |(.Bi+Ci-Ai-.Di)+eisiii0i + e2siii02 

»2 = § (61 + di - 01 -ci + 2a). (17) 

As expected, we see that differential gain errors gi — g 2 lead 
to instrumental Q polarization, 71. 'Spin- flip' errors, cou- 
pling P to P*, are absent at first order but appear at sec- 
ond order in the perturbation and higher. It should be noted 
that the validity of the perturbative expansion depends in 
part on the relative amplitudes of the polarization and total- 
intensity fields. For example, we are implicitly assuming that 
any parameter that contributes to a at first order and to 71 
at only second order is sufficiently small to suppress the 
total-intensity leakage caused to well below the level of the 
polarization leakage. 

For the half-wave-plate receiver, with the plate rotating 
at an angular velocity ip, the ideal Jones matrix is 

Jrhwp — Joint J rot Jhwp J rot 

1 \ / cos ipt — sin ipt 
1 J \ sin ipt cos (fit 



1 

-1 



cos ipt sin ipt 
— sin ipt cos ipt 



cos 2ipt sin 2<pt 
sin 2ipt — cos 2ipt 



(18) 



This leads to similar ideal detector outputs as the pseudo- 
correlator, but with Q and U modulated at a frequency of 
4tp : 

Di = ± (T + Q cos 4ipt + U sin 4ipt) 

D 2 = \(T -Q cos 4ipt-U sin 4ipt). (19) 

Systematic errors in the OMT are parametrized as in equa- 
tion (16), and for the other components, 



Jhwp 



1 + hi Cie lxi 
C 2 e iX2 -(l + h 2 )e lf > 

cos (ipt + 8<p) sin (ipt + 5(f>) 
— sin (ipt + 54>) cos (ipt + 8<j>) 



(20) 



Propagating these errors through to the observed Stokes pa- 
rameters, we find the only non-zero polarization couplings 
are 

a = 9i + 92 + hi + h 2 
2u) — ei cos 61 — e 2 cos 9 2 — 45(f> 

-Cicosxi -C2COSX2. (21) 

The observation that P couples only to P actually holds 
exactly for this receiver, and not just to first order. By com- 
paring equations (17) and (21) we can begin to draw some 



useful conclusions as to the relative suitability of these re- 
ceivers for CMB polarimetry. The half-wave-plate receiver 
has the potentially significant advantage of having no total 
intensity leakage, given the assumptions made. The large 
difference in the amplitude of the temperature and polar- 
ization signals means that such leakage is potentially very 
damaging, and hence any systematic errors that contribute 
to 71 and 72 will have very strict tolerance limits. It should 
be noted that, in the presence of realistic optics, it is likely 
that systematic errors in the half-wave-plate receiver will 
contribute to such leakage, but not at first order, as seen 
for the pseudo-correlation receiver. Also, for a quasi-optical 
wave plate, further investigation is still required into poten- 
tial effects such as optical distortions introduced by the plate 
(Johnson et al. 2007). 



4 BEAM MULLER FIELDS 

The Miiller matrices of the previous section describe prop- 
agation of the Stokes parameters through the receiver. The 
propagation through the entire instrument, including the 
optics, is described by matrix-valued fields. These couple 
the polarization fields on the sky to the observed polariza- 
tion that we assign to the nominal pointing direction. If we 
ignore the effects of reflections from different stages of the 
receiver (and the associated standing waves they set up), we 
can multiply the Miiller matrices for the various stages and 
the optics to get the matrix describing transfer through the 
entire system. In this section we consider the Miiller matrix 
associated with the instrument optics (and horn, if present). 

For a dual-polarization system, let the two orthogo- 
nal polarizations have associated far-field radiation patterns 
EA(n) and iSs(n). It is convenient to define the fiducial 
orientation such that the nominal directions on the sky de- 
fined by E a and Eb (i.e. their co-polar parts) are along the 
x and negative j/-axis of the sky coordinate system, which, 
recall, is the basis we are using to define Stokes parame- 
ters. Of course, for a non-ideal system there will also be 
cross-polar components of Ea and Eb perpendicular to the 
co-polar components. If we adopt the Ludwig-III standard 
for the co- and cross-polar basis, the cross-polar direction 
for A is along — y and for B is along —x. We shall express 
the components of the radiation patterns in terms of co- and 
cross-polar (complex) amplitudes. The coupling of polariza- 
tion A to the incident radiation along a line of sight n is 
oc Ea(ti) ■ E(h) and similarly for B. The coherency of these 
signals defines Stokes parameters that are weighted integrals 
of the Stokes parameters on the sky and from which we can 
extract the beam Miiller fields: 



M T t = 


\(\Ea\ 2 + \E b \ 2 ) 




M T q = 


\(\Eaco\ 2 — \Ea cross 2 + -Es cross 2 — \E B co 


I 2 ) 


Mtu = 


\ (Ea coEa cross ~ Eb co Eb cross) + c.c. 




M T v = 


\i(EAcoE\ cross + Eb co E* B cross ) + C.C. 




M QT = 


\(\E A f- \E B f) 




Mqq = 


\(\Ea<zo\ — \Ea cross + -Ebco — J Eb cross 


I 2 ) 


M QU = 


\ (Ea coEa cross + Eb coEb cross) + C.C. 




M QV = 


\i(E 'a co-E 'Across — ^Bco-Bscross) + C.C. 




M UT = 


i( — EacoE b CIOSS + Ea cross Eg co ) + C.C. 
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— 5 ( — -Ea co-Eb cross — Ea cross Eg co ) + C.C. 
= t;(EacoE b co — Ea cross -Bb cross ) + C.C. 
= \l(E,A co -Es co + -Ba cross -Bfl cross) + C.C. 
= \i{E>A coEg cross — EactossEb co ) + C.C. 
^i(E A co E B CIOSS + i?A cross -Efl co ) + C.C. 



M[7Q 

= jH-^Aco^Bcross "I" -C' A cross -&B co J 
MyU = 2*( _ EacoE Bco + i?A cross -Bb cross) + C -C. 

co Eg co + #A cross S cross 

)+c.c. (22) 

Expressions for the complex Miiller matrices, obtained from 
linear combinations of the components above, are given in 
Appendix A. 

With suitable normalization, the (complex) Stokes vec- 
tor for the signal at the input to the receiver when the beam 
is translated by x is 



p obs (x) = J d 2 x'M{x')p(x + x') 



(23) 



If we first rotate the beam by ip (from x towards — y) and 
then translate by x, the observed signal after back-rotating 
to the sky basis is 

p obs (x;i,) = j & 2 x' A\tP)M{R- 1 x')K{tP) P {x + x'), (24) 

where R^, generates a rotation through tp. In general, the 
behaviour of the observed polarization under rotations of the 
instrument depends not only on the spin of the appropriate 
Miiller matrix terms but also on the beam shapes. 

For theoretical work it is convenient to decompose the 
beam Miiller fields into components that transform irre- 
ducibly under rotation of the instrument. For co-polar main 
beams that are approximately Gaussian, the following ex- 
pansion is particularly useful: 



m(x) = ^^2M mnC j m+n {d x + id y ) m {d x -id y y 



X e 



(25) 



where the sum is over integers m, n ^ 0, and a is the nomi- 
nal beam width. We show in Appendix B that this is related 
to a Gauss-Laguerre expansion and detail the inversion to 
obtain the matrix- valued coefficients M mn . Note that the 
symmetries in equation (8) hold for the irreducible compo- 
nents if we interchange m and n, for example [M mn ]p*T = 
[M nm ] PT . The m — n — components describe pure Gaus- 
sian Miiller fields; their effect is the same as a receiver Miiller 
matrix acting after ideal Gaussian optics. Under a rotation 
of the instrument, we have 

mn 

x A t (i/))M mn A(i/>) 

x (d x + id y ) m (d x — id y )"]e~ x /2<T . 

(26) 

Inserting the beam expansion into equation (24) and inte- 
grating by parts we can express the observed polarization 
in terms of a local expansion in derivatives of the Gaussian- 
smoothed polarization field, p(x; a): 



Pobs(a:;^) 



A t (V)M mn A(i/ ) )(-a) r 



X (d x + id y ) m (d x -id y ) n ]p(x;a). (27) 

The derivative terms have spin m — n plus the intrinsic spin 
of the polarization field, and the coefficients M mn have spin 
n — m plus the difference of spin between the fields involved, 
e.g. the Miiller element [M mn ]pT has spin n — m + 2. This 
expansion generalizes the local expansion introduced by Hu 
et al. (2003). For optical systems where the non-ideal be- 
haviour can be parametrized by only a few M mn matrices, 
equation (27) leads to a fast way to simulate the signal com- 
ponent of polarization maps contaminated by systematics 
for an arbitrary scan strategy (encoded in the angles ip). We 
discuss this further in Section 5. 



4.1 Polarization calibration with beam rotation 

Instrument rotation is a powerful way of reducing the im- 
pact of imperfections in the polarimeter. Structure in M(x) 
that is not invariant under rotation produces systematic ef- 
fects in the maps of observed Stokes parameters that can 
be suppressed by a carefully designed scan strategy, involv- 
ing multiple visits to sky pixels in a range of orientations. 
However, there remain a class of systematics that are not 
averaged out in this way: those for which 



A t (^)M(R- 1 a;)A(V') = M{x) 



(28) 



for all ■0. It is only these contributions to the Miiller matrix 
that would survive in the ideal case of a scan strategy where 
every pixel is visited in all orientations, and for which the 
effective Miiller matrix is 



,(x) = 



/^At W M(R- : 



)A(V)- 



(29) 



If we focus on the element of M that couples a spin-s' field 
into a spin-s one, equation (28) is satisfied by those irre- 
ducible components for which s — s' + n — m = 0, i.e. those 
which are spin-0. Alternatively, if we think of expressing 
the Miiller fields in terms of polar coordinates \x\, 4>, where 
x + iy — \x\e 



[Ms 



,(x) 



2tt 



s)i/> 



M]„'{\x\,4> + i/>) 



i(s-s')4, 



d<fi i(s' 
— e 
2tt 



(30) 

which is just the s — s' (angular) Fourier component. 

As an example, consider the conversion of temperature 
to polarization which is potentially very troubling for CMB 
polarimetry. It is the quadrupole part of M PT (oc e~ 2l< ^) 
that generates instrument polarization that transforms like 
a true polarization under rotation of the instrument, and for 
this we have, using equation (27) 

A^ bs (x) = ' 1 



(d x - idy) 

x ^{Mr, 



1+ 2]pTa 2(m+1) V 2m 2>;a).(31) 



If the summation is real, this is only E-mode contamina- 
tion which can be tolerated at a relatively higher level than 
B-mode contamination. If there is no cross-polarization, 
inspection of equation (A3) of Appendix A shows that 
Mpt(x) is real, and equation (B2) of Appendix B shows 
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that [M m m+2] will then be real provided that the quadrupo- 
lar part of M PT has its planes of symmetry aligned with the 
x- and t/-axes. Under these conditions, the leakage of tem- 
perature to polarization is pure i?-mode (Hu et al. 2003). 

Axisymmetric optical systems have the property that 
their leakage described by Mpr is a pure quadrupole and 
so is not suppressed by rotation (Carretti et al. 2004). To 
see this, we note that axisymmetry and reflection symmetry 
demand that the beam fields generated by a constant field 
h across the beam-defining element (e.g. a horn) are related 
to h by a tensor-valued field that can only be constructed 
from Sij, XiXj and scalar functions of \x\. This implies that 
25 a (n) and _Es(n) can be derived from two radial functions, 
£i(|a:|) and fa ( | as | ) as 



E A = 
Eb = 



£l + |f 2 COS 2(f) 

if 2 sin 2(f) 

— 2 sin 2(f> 
— f l + jf 2 cos 2(f) 



(32) 



For an alternative derivation in the CMB context, see Bunn 
(2006). The PT Miiller matrix element evaluates to 

M PT (x) = R(£i£ a >- ai *, (33) 

which is a quadrupole. Although the leakage cannot be mit- 
igated by instrument rotation, it is pure electric (Carretti 
et al. 2004) since [M m m+2\PT is purely real. Note also that 
the coupling is only to anisotropics in T so a uniform unpo- 
larized brightness would produce no leakage. 



4.2 Optical errors 

In this section we illustrate the ideas developed above by 
analyzing some simple, but common, beam patterns. Taking 
the ideal radiation fields as co-polar Gaussians of width a, 
we consider the true co-polar beams to have pointing errors, 
and ellipticities aligned with the instrument axes, i.e. 



Eac 



1 



1 / (*-l>A,s 
(l + e A ) 

1 ( (*-l>B,x 
(l + e B ) 



s,y> 2 \ 



34) 



These are normalized so their squares integrate to unity. The 
errors can be conveniently reparametrized in terms of the 
average and differential ellipticity and pointing errors (Hu 
et al. 2003), 



P 



b d 



q = 



b B ) 



-(e-A + es) 
i(e A -en). 



(35) 



Cross-polar beam patterns tend to be design-specific and 
are not easily generalized. We gave one example in Sec- 
tion 4.1 where we discussed axisymmetric systems; here we 
consider two simple toy-models instead based on low-order 
quasi-optical approximations. The first model has the cross- 
polar beams as co-pointing Gaussians with the same width 



as the ideal co-polar beam 

VA 



E A cross 
Eg cross 



V2na 2 

VB 
V27TCT 2 



* 2 +y 2 , ■ 



(36) 



where the parameters va and vb control the amplitudes, 
and \a and xs the phases, relative to the co-polar beams. 
In our second model the cross-polar beams have a line of 
symmetry along one axis: 



Eac 



vaV 



V2 



J Bcross — 



*+y A 



V2 



(37) 



In all cases the fields are normalized so the integrals of their 
absolute squares are v A and v B - As with the co-polar case, 
it is useful to reparametrize in terms of average and differ- 
ential quantities, in this case the real and imaginary parts 
of i/Ae' XA and v>Be lXB - the average and differential com- 
ponents of cross polarization in phase and iv/2 out of phase 
with the main beams: 

Va,R = \{~VA COSXA + VB COS^b) 
Vd,R = |(^COSXA - ^SCOSXfl) 

v a ,i = \{va sinxA + v B sinxs) 
Vd,i = \{yAsm.XA — vbsxo-Xb)- 



(38) 



It is cumbersome to proceed exactly in the presence 
of the pointing errors, so instead we expand in the small 
systematic parameters. For the pointing and ellipticity pa- 
rameters, such a low-order expansion is only useful above 
the beam scale a. Even small pointing errors (relative to a) 
and ellipticities can produce large effects below the beam 
scale. This is a clear driver, quite apart from the secondary 
scientific benefits, for increasing the resolution of CMB po- 
larization experiments. The next-generation of ground-based 
B-mode experiments have planned resolution < lOarcmin, 
so a perturbative treatment should be sufficiently accurate 
on the scale of a primordial B-mode signal. 

The Miiller fields follow from equation (22) and their 
irreducible components, M mn , are given in Appendix B to 
first-order in small parameters. (The second-order terms 
that are needed for a consistent polarization power spec- 
trum calculation are also given; see Section 5.) From these, 
we can construct the observed Stokes fields when the instru- 
ment is observing at angle V from equation (27). The result 
for the linear polarization for Gaussian cross-polar beams is 



obs 



= [1 - 2iv a ,R, + \e % ^G\p{d x - id y ) 

+ ^a- lP (d x + idy) + \e- 2 ^G 2 e 3 {d x - id y f 
+ ±e 2 ^o 2 e 3 {d x + idyf}P{o) 



+ [2ie ^vd^p, + 2<s ? <ribd(d x - id y ) 

i 

2 



+ ie 3 "^ 'a -ib d (d x + id y ) + \a 1 q{d x - id y ) 2 
+ ^a 2 q(d x + tdy) 2 ]T(a) 



+ 2e 2i *v d ,iV(a), 



(39) 



where ±\b d = bd ■ {x ± iy) are the spin-±l components of 
b d , and similarly for p. In the fiducial orientation, tp = 0, 
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this reduces to 

Pobs = [l-2w a , R + ap- V + a 2 e s (dl~dy)]P(a) 
+ \2iv d>R + ab d -V + <r 2 q{dl - d 2 )]T(a) 
+ 2v dJ V(o), (40) 

in agreement with Hu et al. (2003) if we drop the cross-polar 
terms. We see that a pointing offset in both beams, p, cou- 
ples to the polarization gradient, whilst a differential point- 
ing error, b d , couples to the temperature gradient. Average 
and differential ellipticity errors, e 3 and q, couple to local 
quadrupole-like patterns in the polarization and tempera- 
ture respectively. As the cross-polar beams have the same 
shape as the ideal co-polar ones, the cross-polar errors cou- 
ple to the Stokes fields directly, and not through gradients. 
Hence, these errors are similar in form to those introduced 
by the receiver, and their effect on polarization can be repre- 
sented by a rotation error, lo = — v a ,R, temperature leakage, 
72 = 2v d , R , and circular polarization leakage, wi = 2v d j. 

For our second, parity-odd, model of cross-polar beams, 
the v terms in equation (39) should be replaced by 

APobs = \e ov at B.{d x - id y ) - e % *crv a , R {d x + id y )]P(a) 
- [e v ov d>R {d x - id y ) - e izv ov d>R {d x + id y )]T(a) 
+ i\e %i 'av d j(d x — id y ) — e i%i ' ov d>I (d x + id y )]V(a), 

(41) 

so that, for ip = 0, the odd-parity cross-polar beams con- 
tribute errors 

APobs = -2iv a , R ad y P(a) + 2i(jv d>R d y T(o) 

+ 2av d jd y V{a). (42) 

That is, the average and differential components of cross po- 
larization in phase with the main beams couple to the gradi- 
ent of linear polarization and total intensity in the direction 
parallel to the line of symmetry, respectively, and the dif- 
ferential component 7r/2 out of phase with the main beams 
couples to the circular polarization in a similar manner. 

For an ideal scan, all terms with ip dependence in equa- 
tions (39) and (41) average to zero. The only terms to remain 
are then 

Pobs = P{ff) - 2ip a , R P{o) + \qa 2 (d x - id y ) 2 T(a), (43) 

where the second term on the right is only relevant for the 
Gaussian cross-polar case. 



5 POWER SPECTRUM ANALYSIS FOR 
TIME-INVARIANT SYSTEMATICS 

In the previous sections we have described the impact of 
various systematic errors on the observed Stokes parameters. 
To assess fully the cosmological consequences of these errors, 
we should propagate them through to power spectra and, in- 
deed, to cosmological parameters. In this section we consider 
the impact of systematics on the P-mode power spectrum. 
Characterizing this spectrum is a major goal of future CMB 
polarization experiments and controlling systematic effects 
will be of critical importance. The case of systematic effects 
whose projection onto Stokes maps can be described as a 
statistically isotropic and homogeneous random process was 



considered in Hu et al. (2003). This is useful for benchmark- 
ing and is probably a reasonable approximation for some er- 
rors such as the pointing jitter/nutation felt by many space 
and balloon-borne experiments. However, the statistical de- 
scription is unrealistic for many other systematic effects. For 
example, optical imperfections can reasonably be expected 
to be time-invariant and their projection onto the map is 
systematic and determined only by the scan strategy. Here, 
we shall concentrate on such time-invariant systematics, and 
develop semi-analytic methods for predicting their effect on 
the power spectra for arbitrary scan strategies. This is use- 
ful for assessing the impact of uncalibrated constant errors 
and to inform methods for removing the effects of calibrated 
errors to a reasonable level. We also consider the raster scan 
(where the instrument is always in it fiducial orientation) 
and an ideal scan, as defined earlier, as special cases for 
which simple analytic results can be found. 

The raster and ideal scan treat each pixel identically, 
and we can define effective Miiller matrix fields that are 
independent of the position on the sky being observed. For 
the raster scan these are just the Miiller fields in the fiducial 
orientation; for the ideal scan they are M S y m . The observed 
Stokes maps are simply the convolution of the Miiller fields 
with the Stokes fields on the sky and so in Fourier space: 

p ohs (l) = 2 7 rM eff (-i)p(0- (44) 

Here, M(i) is the Fourier transform of the Miiller matrix- 
valued field. In terms of the irreducible components, M mn , 
we have 

2ttM(0 = ^(iZa) m+n M m „ e i(m ^ i e- ;2<T2/2 , (45) 

mn 

and for an ideal instrument 27rM(i) = e~' a ^ 2 \ where I 
is the identity matrix. For these simple scans the observed 
fields in the flat-sky limit are realizations of a statistically- 
homogeneous but generally anisotropic process, where the 
anisotropics arise from the contamination fields. Equa- 
tion (44) allows us to calculate the Fourier transform of the 
P-mode map, B(l), since 

B(l) = i (e 2 ^ 1 P(l) - e- 2 ^ 1 P* (I)) (46) 

from equation (2). 

To estimate the power spectrum we use a simple pseudo- 
Ci estimator approach, giving an estimate of the smoothed 
(band-)power spectrum 

= -. [ d 2 l | Pobs (01 2 (47) 

2f aky J b ldlJ b 

where / s k y is the fraction of the sky that has been ob- 
served and b labels the band. We denote the denominator, 
2/sky J, I dl by Nt since it is a measure of the number of 
independent Fourier modes in the band given the sky cov- 
erage. The estimator ignores several important effects, such 
as leakage from E- to P-modes as a result of incomplete sky 
coverage (Lewis, Challinor & Turok 2002), but here we can 
ignore these as their interaction with systematic effects will 
be of secondary importance. (In our simulations below we 
adopt periodic boundary conditions to avoid the complica- 
tions of P-P mixing.) We also assume implicitly that any 
noise bias is removed from the estimated spectrum, and that 
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the effects of smoothing with an ideal, symmetric Gaussian 
beam are taken account of by multiplying by e . 

By averaging over CMB realizations (denoted by angled 
brackets) we find the mean recovered B-mode spectrum, 
(Q^obs)i as a function of the true band-power spectra, Cjf , 
C E , C B and Cj E (all other spectra are taken to be zero as 
the CMB is not expected to be circularly polarized, and we 
assume parity is not violated) . Making use of equations (44) 
and (46), and ignoring the variation of M c g (I) over the radial 
extent of the band, we find 
3 . 

{Cf obs ) = |E / d0 ! (-iy +/ [At(0 i )M eff (-Z)A(0 ( ) 

x F i A t (0 i )M: ff (-OA(0 i )] J y, (48) 

where the sum is over the P and P* elements. We have 
introduced the matrix of true power spectra 

-CP Cf + Cf C?-C B , (49) 
V -Cl E Cf-Cf C E + C B ) 

and have dropped the V Stokes parameter from Stokes vec- 
tors and Miiller matrices. 

For a raster scan, the receiver errors introduced in equa- 
tion (9) give a recovered power spectrum 

{Cb, ob s) = [(l + a) 2 + £(/? + /|)]Cf 

+ [4o; 2 + i(/ 1 2 + /|)]C i f + i( 7l 2 +72)C 6 T . (50) 

Due to the zeroth-order C B term, the perturbation is first 
order in a, but second order in the remaining parameters. 
Therefore, we must be careful with any physical systematics 
that contribute to a only at second order but to u, f or 7 at 
first order, as their resulting power spectrum contributions 
will be of the same order. For the errors due to the optics we 
proceed from equations (45) and (48) and the Muller-matrix 
components given in Appendix B (or directly from the per- 
turbations in the observed fields given in Section 4.2). To 
get a consistent result to second order in the systematic pa- 
rameters, it is necessary to retain second-order terms in the 
isotropic part of Mpp and the hexadecapole part of Mpp». 
Although these are sub-dominant effects in the map domain, 
they produce second-order corrections to B(l) that are pro- 
portional to the true B-modes, and so produce second-order 
effects in the observed power proportional to Cf . For Gaus- 
sian cross-polar beams the result is 

(Cf obs } = [l-e 2 (a 2 l 2 -a 4 l 4 )-q 2 (2-a 2 l 2 -a 4 l 4 /S) 

- b 2 d (2 + a 2 l 2 )/4 - 2(i£ H + </)]Cf 

+ (b 2 d a 2 l 2 /4 + q 2 a 4 l 4 /8 + 2vl R )Cl 

+ ^IrC E . (51) 

Note that, as the perturbation fields are derivatives of the 
Stokes parameter fields for the co-polar errors, the power 
spectrum perturbations couple to the true power spectra 
via polynomials in al. There is no contribution for the co- 
polar pointing error, p, which is as expected as, for a raster 
scan, this simply leads to a global spatial translation of the 
Stokes parameter fields. Note also that there is no generation 
of B-mode power at leading order from E for the raster 
scan if the optics are pure co-polar. Quite generally, with no 
optical cross-polarization, the coupling of linear polarization 



to linear polarization is diagonal in the instrument basis 
and, moreover, is equal for Q and U to first order in beam 
perturbations. (This follows from inspection of the relevant 
Miiller elements in equation (22).) The first-order coupling 
is through the average of the absolute squares of the two 
co-polar beams and, provided this is parity symmetric, will 
not produce B-modes from E for a raster scan. The form 
of the cross-polar errors in equation (51) follows directly 
from that for the receiver case, equation (50), with u) = 
— Va,R, 72 = 2v d ^ R and the second-order map term for a 
given by a = ~{v 2 tR + v 2 j). In the presence of receiver and 
optical errors, the mean observed power spectrum is the sum 
of equations (50) and (51) (with the zero-order Cjf term 
included only once) plus cross-terms between the receiver 
and beam errors; to second-order the cross-terms contribute 

(AC^ obs ) = (2 72 ^d,fl - 2wiv d j + 8ujv a , R )C B 

+ a 2 l 2 f lQ Cf E /2. (52) 

Note that the cross-terms vanish with the cross-polar beams. 
The terms entering with Cjf and C E follow simply from 
making the replacements 72 1— ► 72 + 2v dtR and u) 1— » u) — v a ,R 
in equation (50). The terms entering with Cjf are from the 
second-order isotropic (m = n) parts of the PP and P* P* 
elements of the combined Miiller matrix for the receiver and 
optics (i.e. the product of their matrices). Finally, the Cj B 
term comes from correlating the B-modes produced from 
real E- modes by spin-flip (/1) errors in the receiver with 
those from differential ellipticity acting on the temperature. 

For our odd-parity toy-model for cross-polar beams, the 
v terms in equation (51) are replaced by 

(AC^ obs ) = -{y\ iR + v 2 a j){2 - a 2 l 2 )C E + a 2 l 2 uj tR Ct 

+ {<T 2 l 2 Va,Rbd,x - 2o 2 l 2 V a ,RV d ,R)Cl E 

+ 2a 2 l 2 vl R Cf , (53) 

correct to second-order in the systematic parameters. Note 
the appearance now of the cross-power C^ E that arises from 
E-B conversion due to u a ,R and T-B conversion from bd or 
Vd,R- If we also include the receiver errors, the only cross 
term to arise now is a 2 l 2 qfiC^ E /2. 

As we have seen in Section 4.1, an ideal scan in which 
each pixel is visited in every orientation suppresses system- 
atic errors that are not invariant under instrument rotation. 
The effective Miiller matrix for this scan is M sym , as defined 
in equation (29). Hence, we can calculate the recovered B- 
mode spectrum in a similar manner as for the raster scan, 
but now only contamination fields with the correct spin will 
contribute. For Gaussian cross-polar beams and receiver er- 
rors we find 

(Cbf obs ) = [(1 + a) 2 + 2y2.v d ,R :- 2w x v dJ + &LUv a ,R 

- e 2 (a 2 l 2 - a 4 l 4 /2) -q 2 {2- a 2 1 2 - a 4 l 4 /S) 
-a 2 l 2 p 2 /2-b 2 i {2 + a 2 l 2 )/A 

- 2{ul R + u 2 a j)]C^ 

+ 4(uj-is a , R ) 2 Ci : . (54) 

The qualitatively new effect here is the appearance of the 
pointing error. A fixed (in the instrument frame) average 
displacement of the pointing centre from its assumed posi- 
tion leads to a symmetric beam distortion on averaging over 
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all orientations. For a small pointing error the tendency is 
to increase the effective beam size and so reduce the power 
spectrum below the beam scale. For the odd-parity cross- 
polar case, the ^-dependent terms in equation (54) are re- 
placed by 



(AC ft s obs ) = -(yln + < 7 )(2 - a 2 l 2 )C? . 



(55) 



These results suggest that, at least to second order in the 
parameters considered, the effects of total-intensity leakage 
can be entirely removed from the recovered B-mode spec- 
trum by an appropriately-designed scan strategy. 

However, there are other sources of constraints on the 
scan. In particular, for the ground-based experiments we 
are immediately concerned with, the strategy needs to be 
designed to best target patches of the sky free from fore- 
grounds, and to aid the removal of atmospheric emission 
from the data, which generally involves constant-elevation 
scans. This will limit the amount of cross linking in each 
pixel that can be achieved (i.e. the spread of tp over which we 
can observe each pixel). Since each pixel is treated differently 
for a realistic scan, the map-domain effects cannot be repre- 
sented by a simple convolution with an effective Miiller ma- 
trix. This raises the possibility that systematic errors that, 
for the raster and ideal scans only produce contamination 
proportional to the true B-modes, may cause leakage of the 
E into B due to the non-trivial geometric pattern of the 
scan. That is, certain errors may lead to more significant bi- 
ases when considered with a realistic scan than for either of 
the extreme cases we have so far considered. Note also that 
the observed fields are neither statistically-homogeneous nor 
isotropic for general scans. 

In order to investigate this, we derive the mean esti- 
mated power spectrum for an arbitrary scan strategy, de- 
fined by a list of values of ip for each pixel, using the Miiller 
matrix decomposition introduced in Section 4. Introducing 
appropriate notation, the pixel at x is visited N x times dur- 
ing the scan, and the angle of the instrument basis relative 
to the fiducial basis on the ith visit is ipi(x), where i runs 
from 1 to N x . Averaging over the observations for each visit, 
the observed Stokes vector is 



(56) 



Re-writing this in component form, and using equation (27) 
we have 



Poba,j{x) = ^{(-ff)™ + X-»+ St - Sj W[M 



mn\ jk 



x (d x + id v ) m (d x - id y ) n } Pk (x; a). (57) 

Here, Sj is the spin associated with the jth element of the 
complex Stokes vector, and 

Rs(x) = j^;J2 e MisM*)}- (58) 
i=i 

This leads to a mean estimated power spectrum of 



8 /sky 



1 f^T\ 



(_iy+i' e i ( s j>- s j)<l>i 



^(-l) m+n [M m „] 3fc [[Mw] 3 V 



f d 2 l' ( 

X / ~2% V ^ m - n + s k- s j ( l ~~ 1 )Rm'-n'+s k ,-s j , (I ~ I ) 
X (i<Tl') m+rl+m ' +n ' e i{m-n-m'+ n '+s k -s kl )^, [ F; ,] fcfc ,^j j 



(59) 



where R s (l) is the Fourier transform of R s (x). For a raster 
scan, R s (l) = 2ir8(l) and equation (59) properly reduces 
to equation (48) after summing over m,n,m' ,n' with equa- 
tion (45). If the optical errors are smooth enough that only 
a few terms M mn in the beam Miiller fields are significant, 
equation (59) is an efficient way to calculate the mean ob- 
served B-mode spectrum from a set of orientations ipi(x) 
avoiding the need for simulations. It is one of the main re- 
sults of this paper and we will make further use of it in 
Section 7.1. 

As well as biasing our power spectrum estimates, sys- 
tematic errors will alter their covariance structure. There 
will generally be an increase in the random error in the esti- 
mates, and, for inhomogeneous scans, additional correlations 
in the power spectra above those due to the survey geometry 
and any inhomogeneities in the instrument noise properties. 
These changes may also impact upon the usefulness of an 
experiment. To assess their extent, we can calculate the co- 
variance of the observed spectrum over realizations, 



C0v(C , ( , jObs , Cy >obs ) 



d 2 ld 2 l'\{B oha (l)B obs (l'))\ 2 



bJV 



(60) 



where we have ignored the non-Gaussianity of the lens- 
induced B-modes (Smith, Hu & Kaplinghat 2004; Smith, 
Challinor & Rocha 2006). For the case of an ideal instru- 
ment, the off-diagonal terms vanish and 



var(C,,, obs ) 



Idl 



(61) 



which is the standard result for cosmic variance. If we 
also include the effects of homogeneous white noise with 
a power spectrum , this expression still holds, with 
Cjf i — ► Cjf + iVjf . For the raster and ideal scans it is quite 
straightforward to determine analytically the variance in the 
presence of the specific systematic errors defined in previous 
sections. To avoid over-cluttering our formulae, we ignore 
any cross terms between different systematic parameters so 
the following expressions are only applicable when all but 
one of the parameters are set to zero. These results are given 
to fourth order in the parameters but with some necessar- 
ily sub-dominant terms neglected: the co-efficients of each 
quadratic power spectrum term are only given to leading 
order in each parameter. For a raster scan, the receiver er- 
rors lead to a change in the variance in the recovered power 
spectrum of 



A[var(C ft s obs )] = _{[4a + 3(/ 2 + / 2 2 )](C i f) 2 

JVf, 

+ [16^ 4 + 3(/ 1 4 + / 2 4 )/8](C b I 
+ 3( 7 i+ 72 4 )(C t T ) 2 /8 
+ [8u 2 + (f 2 + fi)]C b B C b E 
+ (7i 2 + 72 2 )CfC b T }, 



(62) 



mnm'n' kk' 



whilst for the optical errors with Gaussian cross-polar 
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beams, the change in variance is 
1 



A[var(C£, bs )] = ±{[4e 2 (2a 4 l 4 -a 2 l 2 ) 



+ q 2 (a 4 l 4 /2 + 4a 2 l 2 -8) 
-(b% x + bi y )(2 + a 2 l 2 ) 

- 8(^a,fl + 

+ 6(4« + 4/)](Cf) 2 
+ 32< fl (Cf) 2 
+ [3 9 V 8 Z 8 /64 + 12*4^ 
+ 7(bX x + bX y )a 4 l 4 /32](C?) 2 
+ [q 4 (2-4a 2 l 2 +5a 4 l 4 /2 
-a 6 f/2 + 5a 8 l 8 /128) 
+ {bd, x + b% y ){8 - 8a 2 l 2 + 3aV)/64 
+ 16i£ B + 2(4, + z4, H )]Cf Cf 
+ [q 2 a 4 l 4 /2 + 

+ q\2a 2 l 2 - 2a 4 l 4 + a 6 l 6 )C b B C? E } ■ 

(63) 

For the odd-parity cross-polar case, the v dependent terms 
in equation (63) are replaced by 

A[var(C 6 s obs )] = ±{[(4^ + 4^(6- 6a 2 l 2 + 9*Y/4) 

JVf, 

+ (4 fl + 4/)(4a 2 / 2 -8)](Cf) 2 
+ I2v 4 a , R o 4 l 4 {Cff 
+ \\vl R o 4 l 4 l2(Cl) 2 



+ [q 2 a 4 l 4 /2 + (bl x + bl y )a 2 l 2 



+ l(vU + v 4 d , I )(2~2* 2 l 2 +3a 4 l 4 /4) 



+ 8v 2 a , R <j 2 l 2 \C b B C b E 
+ 4vi R o 2 l 2 C b B C b T }. 



(64) 



For an ideal scan, only the receiver error terms in a and uo 
remain, and for the optical errors with Gaussian cross-polar 
beams, the change in the variance is given by 

A[var((% S obs )] = -L{[2e 2 (a 4 l 4 -2a 2 l 2 ) 

JVf, 

+ qVZ 4 /2 + 4a 2 / 2 -8) 
-(b'i x + bl y )(2 + a 2 l 2 ) 

- 2(PI +P 2 y)° 2 l 2 
-%{vl,B+vlMC^) 2 

+ 32< fl (qf) 2 

+ 16< H C7f C b E }. (65) 

For the odd-parity cross-polar case, the v dependent terms 
in equation (65) are replaced by 

A[var((% s obs )] = -l(4 fl + 4 / )(4 ( 7 2 P^8)(C ( f) 2 . (66) 

In Section 7.1 we assess the relative importance of the bias 
and the increase in random error. 

We can also investigate the impact of systematic errors 
on the mean and variance of the estimated spectra through 
simulations. For an arbitrary scan, performing simulations 
by directly smoothing simulated maps with the beam in the 
appropriate orientation for each observation is computation- 



ally intensive. For optical errors that can be parametrized 
with only a few low-order irreducible components, a more ef- 
ficient method is to evaluate the summation in equation (27) 
directly with simulated Stokes maps and their derivatives. 1 
Alternatively, the simulated maps can be convolved with the 
real beams in a small number of orientations, ip, and these 
convolved maps can then be interpolated off to implement 
the required scan. This second method requires that the 
beams do not have rapid angular variations. Both methods 
allow one to test the semi-analytic results presented here 
for the power spectrum, and, for the variance, allow calcula- 
tion for arbitrary scans (which is cumbersome to do analyti- 
cally). Having found the observed maps, the power spectrum 
can be estimated using equation (47) and compared to es- 
timates obtained from the input sky maps. By combining 
the results for a large number of sky simulations to estimate 
the mean and covariance of the recovered spectra, equiva- 
lent results to those presented above for the semi-analytic 
approach are found. Importantly, the interpolation method 
provides a means to check the validity of low-order expan- 
sions in irreducible components. 



6 BIASES IN COSMOLOGICAL 
PARAMETERS 

One of the major goals of CMB experiments is to measure or 
constrain cosmological parameters. In the present context, 
the most important parameter is r, and so it is necessary to 
follow systematic errors through to r. As discussed in Sec- 
tion 1, it is often necessary to control systematic errors to 
much better than the random (instrument noise plus sample 
variance) errors in the power spectrum since, typically, many 
power spectrum estimates are combined into relatively few 
cosmological parameters. It is important to note that the ef- 
fects we are considering properly apply to the residual effect 
of systematics after any attempt to remove them coherently 
during data analysis. Dealing with known receiver errors is 
reasonably straightforward in the time domain - for a recent 
example application to real data see Jones et al. (2006) - but 
dealing with optical errors is considerably more difficult. 

We begin by assuming that no statistical correction is 
made in the power spectrum for systematic effects. The sys- 
tematics will then typically lead to parameter biases that 
we can estimate as follows. A simple estimator for r, equiv- 
alent to the maximum-likelihood estimate for a Gaussian 
likelihood, is 



ct 



(67) 



E b (drC*r/* 2 

where the true C b has been written as a sum from gravita- 
tional waves, rdrC® , and weak gravitational lensing, C^^. 
The variances a 2 are our best approximation to var((7 ( f obs ). 
The bias in r is the average shift in position due to system- 
atics, 



A7* SV s — 



E 6 a r cf(Ac 6 s obs )/a; 



E b (drC b B ) 2 /cr 2 



(68) 



1 The derivative maps can be easily formed in Fourier space. 



© 2007 RAS, MNRAS 000, 1-19 



12 Daniel O'Dea, Anthony Challinor and Bradley R. Johnson 



where (A(7^ obs ) is the bias in the power spectrum, and we 
can ignore any contribution from Avar(C7^ obs ) to <r 2 at lead- 
ing order. A large bias in r compared to the ideal random 
error, 



-1/2 



(69) 



would demand that a statistical correction be made in the 
power spectrum for the systematics. In principle, this can 
be done for known (but unremoved) systematic effects if the 
statistics of the true fields are known a priori, either by sim- 
ulation or with the methods developed in this paper. For 
systematic parameters whose values are uncertain, it may 
be possible to fit their effect in the power spectrum with 
a few parameters that can be marginalised over when esti- 
mating cosmological parameters. An example of the latter 
is the treatment of beam uncertainties in the recent WMAP 
analysis (Jarosik et al. 2006). 

However, even if such statistical corrections can be 
made accurately, systematic effects will generally still in- 
crease the random errors in parameters. We can estimate 
this increase by looking at the variance of the estimator in 
equation (67) in the presence of systematics. Linearising in 
obs ), we have 



Aoy 



£ 6(/ drC? d r C$A C Ov(C£ ohs , C£ ob .)/(<7?<7?,) 



2 [E b (drC») 2 /<?Z] 



3/2 



(70) 



The impact of systematic effects will be negligible if both 
Ar sys and Aoy are much less than the 'ideal' random error 



7 APPLICATIONS 
7.1 Finding limits 

In this section we set tolerances on the systematic parame- 
ters introduced in Sections 3 and 4.2 such that the bias and 
increase in random error on r are below some small fraction 
of the ideal random error. We work in the context of next- 
generation ground-based CMB polarimeters, and take the 
fractional threshold to be 10 per cent. This should ensure 
that the sensitivity of an experiment to r is not significantly 
compromised. The implications of these criteria depend on 
the true value of r and also the noise and scan properties of 
the instrument. Here we take r — 0.01, a realistic limit for 
the next-generation of CMB polarization experiments. Of 
course, if the true r is greater than 0.01 the impact of sys- 
tematic effects will be less. Given these criteria, it is useful 
to define, 



Ar s 



and 



(3- 



A<7 r 



(71) 



(72) 



so that we demand \a\ < 0.1 and /3 < 0.1. 

For the properties of the instrument and survey, we 
use values appropriate to upcoming ground-based experi- 



ments like QUIET 2 and Clover 3 . Specifically, we assume 
around 500 useable detectors each with 200-^Ks 1/2 NET, 
and a total integration time of one year over a survey re- 
gion ~ 600 deg 2 . For simplicity, we analyse a single square 
patch of side 25° and adopt periodic boundary conditions to 
avoid geometric mixing of E- and B-modes. We note that 
planned ground-based experiments that are not situated at 
polar sites will necessarily require a less contiguous survey 
geometry than that assumed here, but this will have little 
effect on our conclusions. We take the ideal Gaussian beams 
to have full- width at half-maximum of 10 arcmin and use 
512 2 pixels. The noise is assumed white and the entire array 
gives a one-year sensitivity of 3.3 pK-arcmin. We assume 
that every pair of detectors suffers exactly the same system- 
atic effect which can be thought of as the worst-case limit. 
We adopt a simple cosmological model consistent with the 
recent WMAP three- year analysis (Spergel et al. 2006). Un- 
der these circumstances, if r = 0.01 then ay = 0.005 (with 
no attempt to clean out the lens-induced B-modes) and an 
ideal instrument might expect a 2a 'detection' of such an 
r. We further assume that only primordial fields are present 
on the sky. Our simulations include lens-induced B-modes 
by direct remapping of Gaussian realizations of the primary 
CMB polarization with a Gaussian lensing-deflection field. 

We shall consider three different scan strategies: a raster 
scan, an ideal scan and a semi-realistic scan. The raster scan 
is a useful fiducial scan, and the ideal scan represents the 
best-case scenario. The semi-realistic scan allows us to as- 
sess the likely impact of beam rotation on controlling con- 
tamination fields. Since the construction of a realistic scan 
is very instrument- and site-specific, here we consider only 
a toy-model scan that provides for some variation in the 
angles tp between pixels across the field. The scan path is 
shown in Fig. 2; it is based on a constant-elevation scan from 
Dome-C in Antarctica - a leading candidate to site future 
ground-based CMB experiments. 4 In using this scan, the ef- 
fects of differential sky coverage on the noise are ignored, i.e. 
we assume that each pixel is observed to the same depth. 

We use a combination of the analytic results in Section 5 
and simulations to find the mean change in the power spec- 
trum and its variance and hence find tolerance limits for each 
systematic parameter considered in Sections 3 and 4.2. We 
vary each parameter individually and so do not account for 
cross-terms between parameters. The results are summer- 
ized in Table 1 for our toy-model of odd-parity cross-polar 
optics. The corresponding power spectrum perturbations for 
the semi-realistic scan are shown in Figs. 3-5. For all param- 
eters apart from e s and p, the power spectrum perturbations 
and tolerance limits were evaluated using the semi-analytic 
and both of the simulation methods outlined in Section 5, 
with the different methods showing good agreement. Fig- 
ure 3 compares the power spectrum perturbations obtained 
using the analytic and interpolation simulation method for 
a selection of the parameters for the semi-realistic scan. The 



2 http://quiet.uchicago.edu/ 

3 http://www-astro.physics.ox.ac.uk/research/expcosmology/ 
groupclovcr.html 

4 The Clover experiment was originally planned to be sited at 
Domc-C but this has since changed to Atacama, Chile. The 
BRAIN pathfinder experiment was deployed at DOME-C for the 
2005-06 austral summer. 
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Figure 2. The path of the central focal-plane pixel for the semi- 
realistic scan. The scan is based on a 24-h constant-elevation scan 
from Dome-C. The green (solid) square shows the region used in 
the analysis and over which the fields were taken to be periodic. 



Table 1. Tolerance limits for systematic errors. Three different 
scan strategies are considered: a raster scan, an ideal scan, and 
a semi-realistic scan (see text for details). The cross-polar results 
are for our odd-parity toy-model. 




Figure 3. Power spectrum perturbations for bd, x and q for the 
semi-realistic scan calculated with our semi-analytical method 
(black, solid and blue, dashed respectively) and with 1000 simula- 
tions (orange, dash-dotted and green, long-dashed respectively). 
The error bars are estimates of the theoretical error in the mean 
of the 1000 simulations due to their finite number. Note the good 
agreement between the two methods. For the simulations the scan 
is interpolated at intervals of 3.6°, and it is this interpolation that 
leads to the small deviations between the two methods evident 
for q at I < 200. The B-mode power spectrum from gravitational 
waves with r = 0.01 is shown in red (dotted) for comparison. 



Parameter 


Raster 




Ideal 




Semi-realistic 


a 


5.10 X 10" 


-3 


5.10 X 10 


-3 


5.10 X 10 -3 


UJ 


4.30 x 10' 


-3 


4.30 x 10 


-3 


4.30 x 10" 3 


h 


1.27 x 10" 


-2 


CO 




5.33 x 10" 2 


h 


1.18 x 10" 


-2 


CO 




5.23 x 10" 2 


7i 


1.61 x 10" 


-4 


CO 




2.22 x 10" 4 


72 


1.51 x 10" 


-4 


CO 




2.10 x 10" 4 


e s 


5.6 x 10" 


1 


5.5 x 10" 


1 


3.1 x 10" 1 


g 


1.20 x 10" 


-2 


7.26 x 10 


-2 


5.39 x 10~ 2 


Px 






7.3 x 10" 


1 


1.0 x lO" 1 


Py 






7.3 x 10" 


1 


6.0 x 10" 1 


bd,x 


2.05 x 10" 


-3 


1.15 x 10 


-1 


2.82 x 10" 3 


bd,y 


2.05 x 10" 


-3 


1.15 x 10 


-1 


2.84 x 10" 3 


»a,R 


2.55 x 10" 


-2 


7.26 x 10 


-2 


2.77 x 10~ 2 


Vd,R 


1.02 x 10" 


-3 


CO 




1.41 x 10~ 3 


v a ,l 


7.26 x 10" 


-2 


7.26 x 10 


-2 


7.25 x 10" 2 


Vd,I 


1.12 x 10" 


-1 


CO 




2.31 x 10" 1 



good agreement of these methods validates the low-order 
expansions used in our analytic work, and is representative 
of all the parameters considered. Power spectrum pertur- 
bations for the remaining parameters for the semi-realistic 
scan, calculated with the semi-analytic method, are shown 
in Fig. 4. For e s and p, the assumption that the parameters 
are small breaks down below the tolerance limit, and so we 
rely on the interpolation simulation method to obtain lim- 
its and power spectrum perturbations (shown in Fig. 5) for 
these parameters. 

In every case, the limiting criterion is the bias in r which 
is always at least an order of magnitude greater than the 
increase in the random error, oy. This is a result of the lens- 



ing induced B-modes dominating the variance in the power 
spectrum on large scales. To see this it is useful to write a 
and (5 as 



r 2 V- 



b 



b 



Ob 



(Jb 



C B 

6, prim 



Avar(C* b £ 



(73) 



(74) 



where we have made use of equations (68)-(72) and the fact 
that <7y|r=o.oi ~ r/2 by design. C^ prim is the primordial B- 
mode power spectrum due to gravitational waves. Written 
like this we see that a depends on the ratio of the bias in 
the power spectrum to the primordial spectrum, whilst j3 de- 
pends on the ratio of the change in the variance to the ideal 
variance in a similar manner. As a result, the large contribu- 
tion of gravitational lensing to the B-modes suppresses /3 rel- 
ative to a. Taking an u error as an example, from the results 
of Section 5 we find <AC 6 s obs )/C b s prim = 4co 2 C b E /C b B prim 
whilst Avar((7 6 s obs )/cr 2 = 8tJ 2 Cf/Cf to leading order in 
lu. As Cjf > C^ prim at all scales of interest, a > f3. Note 
that this assumes we cannot access the reionization infor- 
mation: for an experiment with sufficient sky coverage, the 
reionization peak, which is expected to dominate the power 
spectrum at very large scales, will be resolved and so the 
increase in random error will be more significant. 

As expected, the limits for a and uj are independent 
of the scan strategy, as they simply rescale the Q ±iU spin 
states and the errors they introduce transform like a true po- 
larization. The limit for a is more strict than might naively 
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Figure 4. Power spectra perturbations for the semi-realistic scan corresponding to the tolerance limits for: (left) a (black, solid), w 
(green, long-dashed), 7 (orange, dash-dotted) and / (blue, dashed); and (right) v a ,n (black, solid), R (green, long-dashed), u a j 
(orange, dash-dotted) and v^j (blue, dashed). The B-mode power spectrum from gravitational waves with r = 0.01 is shown in red 
(dotted) for comparison. 




Figure 5. Simulated power spectra perturbations for e s (left) and p x (right) for 1000 simulations for a raster (green, long-dashed), ideal 
(blue, dashed) and semi-realistic (black, solid) scan. On the left, the perturbation for p y for a semi-realistic scan is also shown (orange, 
dash-dotted). The B-mode power spectrum from gravitational waves with r = 0.01 is shown in red (dotted) for comparison. 



be expected, as the error not only amplifies the primordial 
spectrum, but also leads to imperfect subtraction of the lens- 
ing contribution. For the raster scan, there are no limits for 
an average co-polar pointing error, p, as a global displace- 
ment has no effect on the power spectrum. For the ideal 
scan, 71, 72, /1, fi, Vd,R and u d j are not constrained since 
the errors they produce disappear when averaged over all 
orientations. 

Comparing the semi-realistic scan to the raster scan, 
many of the tolerances are relaxed as a result of the cross- 
linking, but the improvement is relatively modest. In the 
case of the tightest tolerances, those for 71 and 72 which 
leak temperature into linear polarization, the limits are only 
relaxed by around 40 per cent, which suggests that our abil- 



ity to control these errors with beam rotation is limited. 
Although an ideal scan removes these errors completely, a 
scan that covers 90 per cent of basis orientations only relaxes 
the limits by a factor of 10. That is, the cross-linking needs 
to be almost complete before we see significant changes in 
the tolerances, and so demanding such improvement would 
put considerable constraints on the scan strategy. Such con- 
straints are unlikely to be compatible with other constraints 
such as the desire to use constant-elevation scans to avoid 
gradients in the average atmospheric signal. Having the abil- 
ity to rotate the instrument directly about its boresight, 
rather than relying on rotation induced by scanning (and 
sky rotation), will likely prove very useful for efficient sup- 
pression of some systematics in future ground-based surveys. 
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The semi-realistic scan introduces a further new feature 
over the raster and ideal scan: not every pixel is seen at the 
same angles ijj, so there is differential beam rotation across 
the field. This gives rise to new effects, such as transforma- 
tion of i?-modes into observed _B-modes even for an exper- 
iment with no cross-polarization, and identical coupling to 
the sky and propagation through the receiver for the two 
polarizations. Consider an otherwise ideal instrument with 
equal elliptical co-polar beams for the A and B polariza- 
tions. If the range of ip were the same in every pixel, the 
observed Q and U maps would simply be the true maps 
convolved with some effective beam; this does not produce 
any transformation of S-modes into B. Similar comments 
apply to a common pointing error, p, with otherwise ideal 
optics. The effect of differential beam rotation in the errors 
from the average ellipticity, e s , and pointing parameters can 
be seen in Fig. 5, which shows the power spectrum perturba- 
tions for these parameters for all three scan strategies. For 
both parameters, transformation of _E-modes to _B-modes 
leads to an increase in the bias at low I. However, this does 
not lead to a significant tightening of the tolerances. This 
is partly fortuitous, as the E-to-B errors from differential 
beam rotation contribute with opposite sign to the B-to-B 
errors that the parameters contribute for an ideal scan. The 
E-to-B effect may be more significant for a more realistic 
scan, perhaps requiring multiple elevation scans to cover the 
survey region. (Interestingly, demanding better cross- linking 
to control temperature leakage, and suppress low-frequency 
noise, is likely to result in a more complex scan strategy with 
increased differential rotation.) As an extreme example, we 
can examine the effects of a scan in which the basis orien- 
tation for each sky pixel is selected at random, as we might 
expect this to maximize differential rotation. With such a 
scan, the limits become 0.024 for e s and 0.058 for \p\ which 
are still not particularly tight. These results suggests that, 
provided we avoid pathologically-constructed scans, trans- 
formation of i?-modes to B-modes from differential rota- 
tion is unlikely to be significant. However, transformation of 
E-to-B through cross-polarization always warrants careful 
control. 



7.2 Real Beams 

So far we have considered the impact of specific beam non- 
idealities that can be simply parametrized. However, there 
may be significant effects from beam characteristics that are 
unexpected, or not easily parametrized and are not present 
in our simple beam models. Therefore, in this section we 
will extend our scope to examine the effects of real, or sim- 
ulated, beam patterns. For simplicity, we will only consider 
raster and ideal scans. From simulations of the far-field pat- 
terns, Ea and Eb , we can form the beam Miiller fields using 
equation (22). Then, from equation (48), we can calculate 
the bias in the recovered power spectrum, and from that 
the bias in r. We simulated beam patterns for a realistic 
optical setup with a physical-optics code. (Specifically, the 
simulations were done in the context of the Clover experi- 
ment.) Beams were simulated for three different focal plane 
positions: the central pixel and two pixels at the edge of 
the array, displaced from the centre along the optical axes. 
These pixels are expected to be representative of the entire 
array. A re-calibration of the beam centres and axes was per- 
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Figure 7. Perturbations to the B-mode power spectrum for the 
central focal plane pixel (black, solid) and the edge pixels (blue, 
dotted and green, dash-dotted) for a raster scan. (For the edge 
pixels, the perturbations are negative below the kink at I ~ 200.) 
The B-mode power spectrum from gravitational waves with r = 
0.01 is shown in red (dotted) for comparison. 



formed to remove simulation artifacts. The power profiles in 
co- and cross-polarization for the two polarization states of 
the central pixel and one edge pixel are shown in Fig. 6 
(the remaining edge pixel having similar features to the one 
shown). One new feature, not considered in our previous 
analysis, that is immediately evident from these profiles is 
the subsidiary maxima present in the co-polar beams. These 
sidelobe features would not be easy to parametrize, partic- 
ularly for the pixels at the focal plane edge where the effect 
breaks the azimuthal symmetry of the profiles. When evalu- 
ating systematic effects, we assumed the data analysis would 
be performed with Gaussian co-polar beams, with relative 
powers and beam sizes matched to the simulated beams. We 
analysed each pixel separately, but scaled the noise down to 
a value appropriate for a 500-element array. This is equiva- 
lent to assuming that every pixel in the focal plane has the 
same imperfect optical response. We performed our analysis 
with an ideal receiver. 

For the central focal plane pixel a is around 2 per cent, 
irrespective of the scan strategy. Of the pixels considered, 
this is the worst case (with a well below the per cent level 
for the remaining pixels) and is well within the tolerance lim- 
its we have previously suggested. Furthermore, this bias is 
mainly a result of our assumption of simple Gaussian beams 
in the analysis, and could be reduced by a more sophisticated 
treatment of the average, symmetric co-polar beam profiles. 
Fig. 7 shows the perturbations to the power spectrum for 
a raster scan. The perturbations are largely independent of 
the scan strategy, and due almost entirely to the true B- 
modes: the contribution from total intensity and _E-modes 
is always below I0 -6 fj,K 2 . These results were confirmed by 
direct convolution of simulated Stokes maps with the Miiller 
beams. These simulations also confirmed that the increase 
in the random error in r is sub-dominant to the bias, as 
expected from Section 7.1. 
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Figure 6. Simulated (power) beam patterns for a realistic optical setup. Panels (a)-(d) show i?Aco| 2 > \Ebco\ 2 , I A cross 1 2 an d \Eg cross | 2 
for a central, on-axis pixel, and panels (e)-(h) show the same for a pixel at the edge of the focal plane array. The average total power 
in the cross-polar beams, relative to the corisponding co-polar beams is 2.2 X 10 — 7 for the central pixel, and 2.3 X 10~ 6 for the shifted 
pixel. The contours are at —3 dB intervals, down to —57 dB for the co-polar beams, and —30 dB for the cross-polar beams. 
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8 DISCUSSION 

We have presented a framework to describe, in a general 
manner, the scientific impact of instrumental systematic er- 
rors in CMB polarimetry experiments surveying small sky 
areas. A major focus was the behaviour of systematic errors 
under instrument rotation. We introduced spin-weighted 
Miiller matrices to describe the propagation of the Stokes 
parameters through the receiver to simplify the transfor- 
mation properties of the systematic effects under rotations. 
The optical coupling of the receiver to the incoming radi- 
ation field is described by Miiller matrix-valued fields. We 
gave expressions for these in terms of the vector- valued field 
pattern of the antenna, and introduced a convenient decom- 
position of the matrix fields into components that transform 
irreducibly under rotations. The decomposition is useful for 
analytical work and also provides an efficient means of sim- 
ulating the effect in Stokes maps of systematic effects for an 
arbitrary scan strategy. 

We compared two popular receiver architectures for 
CMB polarimetry and showed that modulating the polariza- 
tion with a half-wave plate has the advantage over a pseudo- 
correlation receiver of producing no instrumental polariza- 
tion for ideal optics. It should be noted that, in the presence 
of realistic optics, the half-wave-plate receiver will contribute 
to instrumental polarization, but at a higher order in the pa- 
rameters describing systematics than the pseudo-correlation 
receiver. This result holds for a waveplate in waveguide 
- further analysis is required for a quasi-optical approach 
where the waveplate is before the beam-defining element 
(and typically is at the beam waist). We also analysed sim- 
ple parametrized models for imperfect co- and cross-polar 
beams and presented the irreducible components of their 
Miiller fields. 

We presented an efficient, semi-analytic method for cal- 
culating the bias in the power spectra from time-invariant 
instrument systematics described by an arbitrary Miiller re- 
sponse, and for a general scan strategy. This should prove 
useful for setting tolerance limits during the design phase 
of an experiment. However, since the analysis is entirely in 
the map domain, detailed simulations will still be required 
for data analysis to account properly for the interaction 
of non-trivial map-making with systematic effects. We also 
analysed the extreme cases of scans with no and complete 
beam rotation and obtained simple results for the biases in- 
duced in the B-mode power spectrum, and the increase in 
the random error in the power spectrum. Several important 
effects are missed by these extreme scans, such as the gen- 
eration of B-modes from S-modes, even in the absence of 
cross-polarization, by differential beam rotation across the 
observed field. It is often important to propagate the effects 
of systematic errors through to cosmological parameters to 
assess properly their effect on the scientific returns. We did 
this for the tensor-to-scalar ratio for a survey representa- 
tive of that from next-generation ground-based instruments, 
and set tolerance limits for simple parametrized forms of 
the receiver and optical systematics. The formalism can also 
be applied to non-parametric, simulated or measured beam 
models and we illustrated this with the results of physical 
optics simulations. 

Although instrument rotation is potentially a powerful 
way to mitigate against stable instrumental systematic ef- 



fects, we showed that for a typical scan that can be obtained 
reasonably from the ground, the limited range of angles at 
which a pixel is revisited limits the extent to which system- 
atics can be controlled. For the semi-realistic scan considered 
here we found that tolerance limits on systematics could be 
relaxed by between 40 and 400 per cent over those for a 
simple raster scan. For all cases considered here, the toler- 
ance limits for a realistic scan are never significantly stricter 
than for a raster scan. The distribution of scan angles tp can 
be improved considerably for constant-elevation scans by in- 
cluding z-axis rotation, rather than relying on sky rotation. 
This may prove to be an important element of future cam- 
paigns to detect primary B-mode polarization in the CMB. 

As previously noted, this paper does not directly con- 
sider some important time-invariant systematic errors. For 
example, as we have restricted our analysis to quasi- 
monochromatic systems, errors such as bandpass mismatch, 
and other variations in instrumental response with frequency 
have been ignored. However, the methods and analysis in- 
troduced in this paper can be straightforwardly extended 
to accommodate such errors by splitting the bandwidth up 
into a number of sub-bands in which the instrumental re- 
sponse can reasonably be modelled as constant, and allowing 
systematics to vary between the sub-bands. Each sub-band 
can be treated independently and the Miiller matrix used 
to describe the whole instrument simply becomes an appro- 
priately weighted sum over the matrices for each band, and 
the analysis can proceed unhindered. 
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APPENDIX A: COMPLEX BEAM MULLER 
FIELDS 

We can extract the complex beam Miiller fields from the 
real components in equation (22). Those fields for which the 
second index has non-zero spin are best expressed in terms 
of the spin-±l components of E A and Eb'- 



±iEa = E A ■ (x ±iy) = E Aco =f iE Across 
±\E B = E B ■ (x ± iy) — -E B cross T iE B co, 



(Al) 



whereas those with spin-0 component can be expressed di- 
rectly in terms of scalar products of Ea, Eb and their duals 



*E A 



Eact 

E a , 



*E< 



—Eb co 

Eb cross 



(A2) 



which are obtained by a right-handed rotation through 90° 
about the radiation propagation direction. The independent 
components of the complex Miiller fields are 



Mtt = 


\{\E A \ 2 + \E B \ 2 ) 


M TP = 


j(iEa-iEa + iEb-iE b ) 


Mtv = 


\i(*E A - E* a + *Eb ■ E B ) 


M PT = 


\{E A + iE B )- (E A -iE B y 


M PP = 


\{iE A + iiE B )(-iE A - i-iE B y 


M P p* = 


±(-iE A +i-iE B )(iE A -i 1 E B y 


Mp V = 


\i{*E A + i*E B ) ■ (E A - iE B )* 


M V T = 


\i{-E A - E'b + Eb ■ E A ) 


M VP = 


— ^i(\EA-iE B — \Eb -\E* a ) 


M V v = 


\{*E A ■ E* B — *Eb ■ E* A ). 



(A3) 

Expressed this way, the transformation properties under a 
rotation of the coordinate basis on the sky are manifest. 



APPENDIX B: BEAM EXPANSION 

The expansion of the beam in equation (25) can be written 



M(a;) = 



-a! 2 /2o 



27TCT 2 



■n 



M mn (-l) n 2 m m 



xL" m m (x 2 /2a 2 )], 



(Bl) 



where L™ is the Laguerre polynomial (Gradshteyn & Ryzhik 
2000), familiar from the radial wavefunctions of the Hydro- 
gen atom. Each term in the sum is a Gaussian multiplying a 
polynomial in x and y. Equation (Bl) can be inverted using 
the orthogonality of the Laguerre polynomials: 



±£ J d 2 xM(^Y~ m L" m --(x 2 /2* 2 ).(B2) 



Note that, for m > n with m and n > 0, L r ^ m {x) = 
0(x m ~ n ) so that the integrand in this equation is always 
regular at the origin. 

In Section 4.2 we introduced a simple parametrization 
for offset, elliptical co-polar beams and two toy-model exam- 
ples of cross-polar beams. Here we give the M m „ matrices for 
these parametrized beams to first order in the parameters. 
Generally, if the maximum order of the polynomials appear- 
ing in equation (Bl) is I, then there are terms present with m 
and n taking all integer values > such that m + n < I. For 
a first-order expansion of the beam models of Section 4.2, 
/ = 2. Adopting a matrix notation for the indices m and n, 
we find for the case of Gaussian cross-polar beams that 



-(ip/2) e s /2 



[M mn 


TT = 


"(-1P/2) 












V &s/2 












/ ~Wd,R 




q/4 


[Mmn 


TP = 














V 9/4 
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-(ibd/2) 


q/2 


[M mn 


PT 


-(-ibd/2) 
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/ 1 - 2iv a ,R 


"(ip/2) 


e s /2 


[M m „ 
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-(-iP/2) 












V e s /2 












/ 1 


-(iP/2) 


e 3 /2 


[M mn 


VV = 


-(-iP/2) 












V e s /2 









(B3) 



the remaining independent, non-zero elements arc 
[Moo]tv = 2v a ,R = [M 00 ]vT, [M o]pv = 2u dJ = 
— 2[Moo]vp and [Moo]vt = 2v a ,i- Here, for example, 
±ip = p ■ (x ± y) are the spin-±l components of p. For the 
power spectrum analysis in Section 5 we need some elements 
of Mpp and Mpp» to second-order in the systematics as 
these appear multiplying the zero-order Miiller matrix. The 
relevant second-order, non-zero terms are 



[Mo()]pp 
[Mu] PP 
[M 22 ] pp 
[M 04 ] pp » 



= -&d/4 - q - [v a ,R + v a ,i) 

= p 2 /4 + b 2 d /8 + e 2 /2-q 2 /2 

= e 2 /4 + q 2 /8 

= q 2 /l6= [M 40 ] PP ,. 



(B4) 



© 2007 RAS, MNRAS 000, 1-19 



Systematic errors in CMB polarization measurements 



For the cross-polar beams with a single line of symmetry, 
[M mn ] T T and [M ra „]yy are unchanged from equation (B3) 
but 

1 / ib d + 2v d>R -q \ 

[M mn ] TP = -- -\b d -1vd,R 

V -q / 
1 / i6 d - 2^ d>fl -q \ 

[M mn ] PT = -- -ib d + 2v d , R 

V -q J 

1 / -2 ip + 2f ,/i -e s \ 
[M mn ] pp = -- _ip-2^ a ,H 

2 V -e. / 

(B5) 

The remaining independent non-zero elements are 
[M i]tv = -Wo,/ = -[Mi ]tv, [M i]pv = —Wd,i = 
-[Mio]pv, [M i]vt = -iv a j = -[M 10 ] V t and 
[Moi]vp = ivd,i/2 = — [Mio]vp. The second-order 
terms required for a consistent power spectrum calculation 
are, in this case, 

[M 00 ] PP = -&d/4 - q 2 - {yl tR + v 2 a j) - ip y v a ,R 
[Mn] pp = p 2 /4 + b d /8 + ei/2-q 2 /2 

- (fl,R + V 2 aA)j2 - ipyV a ,R/2 

[M 22 ] pp = e 2 /4 + q 2 /8 
[M 04 ] PP . = q 2 /l6= [M 40 ] PP ». (B6) 
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